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Abstract. As most exoplanets host stars, HD 160691 (alias p Ara) presents a 
metallicity excess in its spectrum compared to stars without detected planets. 
This excess may be primordial, in which case the star would be completely over- 
metallic, or it may be due to accretion in the early phases of planetary formation, 
in which case it would be overmetallic only in its outer layers. As discussed in a 
previous paper (Bazot and Vauclair 2004), seismology can help choosing between 
the two scenarios. This star was observed during eight nights with the spectro¬ 
graph HARPS at La Silla Observatory. Forty three p-modes have been identified 
(Bouchy et al. 2005). In the present paper, we discuss the modelisation of this 
star. We computed stellar models iterated to present the same observable param¬ 
eters (luminosity, effective temperature, outer chemical composition) while the 
internal structure was different according to the two extreme assumptions : origi¬ 
nal overmetallicity or accretion. We show that in any case the seismic constraints 
lead to models in complete agreement with the external parameters deduced from 
spectroscopy and from the Hipparcos parallax {L and Tefi). We discuss the tests 
which may lead to a choice between the two typical scenarios. We show that the 
“small separation” seem to give a better ht for the accretion case than for the 
overmetallic case, but in spite of the very good data the uncertainties are still 
too large to conclude. We discuss the observations which would be needed to go 
further and solve this question. 
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1. Introduction 

Since the first discovery of a planet orbiting around Peg 51 (Mayor et al. 1995), more 
than 130 exoplanets have been detected The central stars of these planetary systems 
appear to be overmetallic compared to the Sun, at least in their atmospheres. Their 
average metallicity is ~ 0.2 dex larger than solar while it is about solar for stars which 
have no detected planets, in the solar neighborhood (Gonzalez et al. 1998, Santos et al. 
2003). Two scenarios have been proposed to explain these high metallicities. In the first 
scenario, they are the result of a high initial metal content in the proto-stellar gas : in a 
metal-rich environment the mechanisms leading to planet formation are more efficient, 
which would explain the observed bias. In the second scenario, the overmetallicity is due 
to the accretion of hydrogen-poor matter during planetary formation. This matter could 
come from already formed planets which migrate to wards the central star (Trilling et 
al. 1998) and eventually fall on it (Murray et al. 2001, Murray & Chaboyer 2002). This 
process leads to a metal excess in the stellar outer layers only (as it is diluted in the 
surface convective zone), while the former scenario implies that planet-hosting stars are 
completely overmetallic, down to their centers. 

Bazot & Vauclair (2004, hereafter BV04) studied the internal differences of stellar 
models computed with either the overmetallic or the accretion assumption, but iterated 
so as to present exactly the same observable parameters (Teg, ■^, Z) : these models 
can account for the same observed star, but their interiors are different. Moreover, they 
have slightly different masses (up to 20%) while their observed parameters are the same. 
BV04 studied the asteroseismic signatures of these internal differences. They focused on 
the fact that something special occurs for stars with masses around 1.1 Mq, as convective 
cores begin to develop. This is a delicate and metallicity-dependent process (Michaud et 
al. 2004). Indeed, it is possible that among two models representing the same observed 
star, one has a convective core while the other does not. 

We decided to test these preliminary investigations with the star /r Ara (HD160691, 
HR6585, GJ691). It is a G3IV-V star with a visual magnitude V=5.1. We obtained 8 
nights in June 2004 to observe this star with the HARPS spectrograph on the 3.6-m 
telescope in La Silla. We first reported on a discovery of a new low-mass planet around 
this star, of ~14 earth masses (Santos et al. 2004b). The star itself presents very clear 
oscillations : the analysis of the Fourier transform of the radial velocity time series allowed 
us to detect 43 oscillation modes of degrees 1=0 to 1=3 (Bouchy et al. 2005, hereafter paper 
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Table 1. Effective temperatures and metal abundances observed for fi Ara. [Fe/H] ratios 
are given in dex. The references for the given values are given in column 3, they are 
followed by f for spectrometric measurements and | for photometric measurements. 


Teff(A) 

[Fe/H] 

Reference 

5798 ± 33 

0.32 ± 0.04 

Santos et al. 2004af 

5813 ± 40 

0.32 ± 0.05 

Santos et al. 2004bf 

5800 ± 100 

0.32 ± 0.10 

Bensby et al. 2003/ 

5811 ± 45 

0.28 ± 0.03 

Laws et al. 2003/ 

5800 ± 90 

0.28 ± 0.04 

Favata et al. 1997/ 

5570 ± 70 

0.39 ± 0.10 

Laird et al. 1985]: 

5597 ± 160 

0.41 ± 0.15 

Perrin et al. 1977/ 


I). A summary of these observations is given below, in section 2. The model computations 
and comparison with the observations are presented in section 3. The results are discussed 
in section 4 while we conclude and propose future work in section 5. 

2. Observational constraints 

With a visual magnitude of V=5.1 and a Hipparcos parallax of 7r=65.5±0.8 mas, the star 
fi Ara has an absolute magnitude My=4.20 and a luminosity of log L/Lq =0.28±0.012. 
The different values given in the literature for its effective temperature are listed in 
Table ^ The HARPS spectrum obtained during our run (see below) was used to de¬ 
rive a new spectroscopic set of atmospheric parameters. From this study the values 
Teff=5813±40 K and [Fe/H]=0.32±0.05 dex were found (second line of Tabled). With 
such a metallicity, /i Ara happens to be one of the most metallic planet-hosting star (see 
e.g. Santos et al. 2003). As we can see in Tabled there is a discrepancy in effective tem¬ 
perature between the two oldest measurements (bottom lines) and the most recent ones 
which are in good agreement, around 5800 K. In Figures ^ and |21 the five most recent 
measurements are displayed. We will see below that seismic studies of this star lead to 
more precise constraints on Tea, inside the error bars given by these spectroscopic values. 
We also notice that they are in good agreement with the newest effective temperatures, 
which are believed to be more precise than the oldest ones. The two lowest values given 
in Tabled for Tea are ruled out by our results. 

We observed n Ara with the high-precision echelle spectrograph HARPS available 
on the 3.6-m ESO telescope at La Silla, Chile. As discussed in details in paper I, 43 
low degree p-modes with a signal-to-noise ratio greater than 2.5 were identified. The 
degrees of the observed modes range from 1=0 to 1=3. The uncertainty on each mode 
frequency is 0.78 ^Hz due to the duration of the time series. The eigenfrequencies for 
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these p-modes are given in Table 2 of Paper I. The average separation between modes of 
consecutive order is Ar'g = 90 /rHz (see paper I). It should be noted that in the further 
analysis we adopt a conservative point of view. Because of the different confidence levels 
we have on the observed frequencies (see paper I for a complete discussion), and thus on 
the individual large separations, we defined the region of seismic constraint as the area 
delimited by the iso-Aizg 89 /rHz and 91 /xHz (see Figs. ^ and 

The echelle diagram for the observed frequencies is shown in both Figures 3 and 4. Its 
interpretation and the comparisons with theoretical models are discussed in the following 
section. 

3. Models 

3.1. Echelle diagrams 

We computed series of “overmetallic” and “accretion” models which could account for 
the observed parameters of g Ara. We used the Toulouse-Geneva stellar evolution code, 
with the OPAL equation of state and opacities (Rogers & Nayfonov 2002, Iglesias & 
Rogers 1996) and the NACRE nuclear reaction rates (Angulo et al. 1999). In all our 
models, microscopic diffusion was included using the Paquette prescription (Paquette et 
al. 1986, Richard et al. 2004). The treatment of convection was done in the framework 
of the mixing length theory and the mixing length parameter was adjusted as in the Sun 
(a = 1.8, Richard 1998, see Richard et al. 2004).The accretion models were computed 
with the same assumptions as in BV04 (instantaneous fall of matter at the beginning of 
the main sequence and instantaneous mixing inside the convection zone). Neither extra¬ 
mixing nor overshoot were taken in account in the present computations (this will be 
done in future work). 

Adiabatic oscillations frequencies were computed for a large number of models along 
the evolutionary tracks, using a recent version of the Brassard et al. (1992) code, as in 
BV04. In order to compare them with the modes observed in g Ara, the frequencies were 
computed for angular degrees 1=0 to 1=3 and radial orders ranging typically from 4 to 
100. The azimuthal order is always m=0. For each evolutionary track, we focused on the 
models which fell inside the observed box in the HR diagram (cf Figures 1 and 2), namely 
models with luminosities and effective temperatures as given by spectroscopy. We then 
searched for those models, among the previous ones, which could precisely reproduce the 
observed echelle diagram. 

For each track, only models with computed large separations of exactly 90 /rHz could 
fit the observed echelle diagram. A difference of only 0.5 p,Hz in the large separations 
would already destroy the concordance. The oblique solid lines in Figures 1 and 2 repre¬ 
sent the location in the HR diagram of the models which correctly fit the echelle diagram 
and lie in the spectroscopic error box. The two other curves (dot-dashed ones) indicate 


Bazot et al.: Asteroseismology of fi Arae 


5 



Fig. 1. This graph represents the constraints on the position of the star fx Ara in the HR 
diagram. Various observational constraints are shown. The error box corresponds to the 
luminosity derived from the Hipparcos parallax and the effective temperature measured 
by Bensby et al.(2003). Other temperature determinations are displayed above the box 
and labeled as follows : 1 - Santos et al. (2004a), 2 - Santos et al.(2004b), 3 - Laws et 
al.(2003) and 4 - Favata et al.(1997). An example of evolutionary track is drawn for the 
“overmetallic” assumption (dashed line, see text). The asteroseismic constraints appear 
as three oblique lines in this HR diagram : the (solid) center line represents the models 
for which the large separations are exactly 90 /rHz as observed in /i Ara, while the two 
other lines stand for 89 /iHz (top dot-dashed line) and 91 /rHz (bottom dot-dashed line). 
The cross shows the place of the model that we have chosen to discuss in more details, 
for the comparison with the accretion case. 


models for which the large separations are respectively 89 and 91 /rHz, which do not fit 
the echelle diagram anymore. 
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Fig. 2. Same as Fig. Q] for the observational constraints on /r Ara. Here the evolutionary 
track is drawn for the “accretion” assumption. The seismic constraints are represented 
for the accretion case in the same way as in Fig. 2] for the overmetallic case. 


We present here two models which, given our basic assumptions, fit reasonably well 
the observed frequencies and echelle diagram. We computed a reduced for both cases. 
The echelle diagrams are given in Figures 3 and 4, respectively for the overmetallic 
and the accretion cases. These models correspond to the crosses in Figures 1 and 2. 
Their characteristics are given in Table 2| The accretion model is labeled AC and the 
overmetallic one OM. The internal distributions of elements in these models are displayed 
in Figured : in the AC model the overmetallicity lies only in the outer layers while it 
is present in the whole OM model. The initial helium mass fraction Yq is larger for the 
overmetallic model than for the accretion model as we assume here that the original 
interstellar matter was overabundant in metals and helium according to the chemical 
evolution of galaxies (Izotov & Thuan 2001, see BV04 for a discussion of this assumption). 
From spectroscopy and photometry both stars would appear nearly identical, within the 
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Fig. 3. Echelle diagram for the overmetallic model. The lines represent the theoretical 
frequencies, the solid line is for 1=0, the dotted line for 1=1, the dashed line for 1=2 
and the long-dashed line for 1=3. The symbols represent their observational counterpart. 
Open squares, full squares, open triangles and full triangles stand for 1=0, 1=1, 1=2, 1=3, 
respectively. Their size is proportionnal to the signal to noise ratio (see Paper 1). 


observational uncertainties, but their internal structures are quite different. The most 
striking fact is the presence of a convective core in the overmetallic model, with a radius 
of ~ 0.07 R* (R* being the radius of the star). 

The accretion scenario, with x2 = 4.3, provides a better match for the observations 
than the overmetallic model, for which we obtained x2 = 7.6. We must however keep in 
mind that these values are still high : it is not possible to give a clear conclusion at the 
present stage. We discuss below other possible seismic tests. 

Among these frequencies, we note that the lowest one is particularly interesting. It has 
been identified as a 1=2 mode (see Paper I). It is in good agreement with the accretion 
model. However, the overmetallic model gives a slightly better fit if we consider it as being 
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Fig. 4. Same as Fig. 01 for the accretion model. Only the mode with the lowest frequency 
differs from Fig. 0 (see text). 

a 1=0 mode. It is a good example of the uncertainty we have on the mode identification 
at very low signal to noise ratio. 

3.2. Small separations 

We compared the measured small separations, Sv = 5 vn,i — r'n-1,1+2, with the computed 
ones. As this quantity is very sensitive to the deep stellar interior (Tassoul 1980, Roxburgh 
& Vorontsov 1994), it was proposed in BV04 as a good indicator of the presence of a 
convective core in stars around 1.1 Mq. 

According to Roxburgh and Vorontsov (1994), the presence of a convective core leads 
to an oscillatory signal in the small separations with a period scaling approximately as 
^ l/ 2 tcc, where tec is the acoustic radius of the convective core (see Table|2l. 

Figures 0 and [7| display the observed and computed small separations for the Z = 0, 2 
and Z = 1,3 cases. We can see an important difference between the overmetallic and the 
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Fig. 5. Internal distribution for elements H, He, C, N, O, Ne and Mg for models OM 
(dashed lines) and AC (full lines). The lower right panel displays the global fraction, Z, 
of elements heavier than helium. 


accretion cases for the I = 0,2 small separations. In the overmetallic case they rapidly 
decrease at large frequencies and even become negative at some point, while in the 
accretion case the curve is much flatter. This behavior has to be related to the echelle 
diagrams : for the overmetallic case (Figure 3), the theoretical lines for Z = 0 and I = 2 
cross around v = 2.7 mHz, which does not happen for the accretion case. Such a crossing 
is in contradiction with the asymptotic theory which may be a good representation of 
seismology in a radiative interior but does not hold in the presence of a convective core. 
It assumes in particular that the sound velocities at the turning points of the waves 
are identical to the central sound velocity, which is not true in this case (Figure IHJ. A 
complete discussion of the signatures of convective cores will be given in a forthcoming 
paper. Here we only compare the observed and predicted small separations in the star /i 
Ara. 
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Table 2. Mass, age, initial and actual chemical compositions, radii, geometrical and 
acoustic, of the star (R* and T*) and of the convective core (Rcc and t^c, when necessary) 
for stellar models which satisfy to the observational constraints for /r Ara. The model 
denoted AC includes accretion, the other one, denoted OM, is overmetallic. Initial mass 
fractions (Xq, Yq, Zq) on the main sequence are given for the core and the envelop, 
which has been enhanced by pollution in the AC case. For the final model, the surface 
abundances are given (X, Y, Z). 


Model 

M* 

Age 

Tefr 

L/Lo 

Xo 

Yo 

Zo 



(Mo) 

(Gyr) 

(K) 


envel. 

core 

envel. core 

envel. 

core 

AC 

1.085 

5.803 

5825.69 

1.887 

0.6913 

0.7097 

0.2650 

0.2714 

0.0437 

0.0189 

OM 

1.180 

4.251 

5769.28 

1.908 

0.6696 

0.6696 

0.2891 

0.2891 

0.0413 

0.0413 


X 

Y 

Z 

R* 

Rcc 

T* 

tec 




(Ro) 

(R*) 

(s) 

(s) 

0.7492 

0.2124 

0.0384 

1.35 

- 

5258.54 

- 

0.6694 

0.2892 

0.0414 

1.38 

0.07 

5240.80 

137.19 


We also computed a reduced x2 for the small separations. We obtained x2 = 1.2 for 
model AC and y2 = 2.2 for model OM. These results are better than those resulting 
from the fits of the frequencies (echelle diagrams). This is partly due to the fact that the 
uncertainties on small separations are times larger than for individual frequencies. 
But there is also a physical reason for this better fit : the small separations are believed 
to be independent of the surface effects which affect the individual frequencies (see e.g. 
Gabriel 1989). These effects, which are difficult to estimate, may deteriorate the quality 
of the fit of individual frequencies while they have no influence on small separations. 

The fact that model AC better matches the observations than model OM has to be 
considered cautiously. We can see in Fig. |H1 that the decrease of the theoretical small 
separations due to the convective core, in model OM, becomes more important at high 
frequencies : this effect, which is crucial for the choice between models, corresponds to 
the beginning of the oscillations that we have discussed above. However, the three values 
that we obtain for the highest observed frequencies involve modes with a signal to noise 
ratio between 2.5 and 3 (small dots in Figs|niand[7I). We cannot exclude misidentifications 
of these modes. In such a case the corresponding small separations would be different 
and the agreement between the observations and model OM could become comparable 
to the one with model AC. 
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Frequency [mHz] 


Fig. 6. Small separations for model OM. The full and open circles are the observed 
values obtained for 1=0-2 and 1=1-3. The small dots represent values involving modes 
detected at a low signal to noise ratio (see Paper I). The theoretical computations are 
represented as a solid curve for 1=0-2 and a dotted curve for 1=1-3. 

More observations are needed to reach a clear conclusion. It would be particularly 
important in this framework to detect modes with large frequencies, even larger than 
in the present work, and with the highest possible signal-to-noise ratio, to be able to 
distinguish clearly between the two scenarios. 

4. Uncertainties 

Stellar models are characterized by their mass (M), an initial helium mass fraction (Yq) 
and a mixing-length parameter (a), and an age. 

Comparing Figs. Hand El it is clear that the region of seismic constraint depends of 
the physics used in the model. The lines of constant average large separation are given 
for a set (a, Yq). Here a is adjusted on the solar models (Richard et al. 2004) and Yq is 
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Frequency [mHz] 


Fig. 7. Same as Fig. |^for model AC. 

derived from the cosmic Y-Z relation (Izotov & Thuan 2001) Changes in a or Yq would 
also lead to modifications of the location of the iso-At'o lines in the HR diagram. To 
estimate this effect, we have computed models with the same input parameters as AC 
or OM except that we imposed a small mixing-length or initial helium mass fraction 
variation (only one at the time). We then selected, for given physics (i.e. accretion or 
overmetallicity and a given set of parameters (a, Yq)), the age for which the average 
large separation is 90 /iHz. 

Our results are given in Table 13 The given variations in Teff and L/L 0 are the 
differences between the values of the recomputed models and those of AC or OM. They 
are representative of the accuracy we can obtain with asteroseismology on a and Yq. 
The changes we imposed on the mixing-length parameter and on the initial helium 
mass fraction correspond respectively to 3% and 4% of the initials values used in AC 
and OM. In fact, variations in Yq with such an amplitude are ruled out because the 
resulting position of the star in the HR diagram would not account for the spectroscopic 
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R/R. 


Fig. 8. Sound speed profiles inside AC (solid line) and OM (dashed line) models. The 
feature in the central region for the OM case compared to the AC case is due to the 
convective core. The sound velocity drops dramatically at its boundary and goes up 
again inside. The frequencies of the modes I = 0 which progress down to the stellar 
center are affected by this behavior while those for which the turning point lie outside 
this region are not. This explains the special behavior of the corresponding lines in the 
echelle diagram also seen in the small separations. 


measurements anymore. On the contrary changing a of 3% has a very small influence 
on the star luminosity and does not affect the effective temperature in such a important 
way as a variation in Yq does. 

The problem is different for the precision obtained on the mass and the age of the 
star. In this case, a variation in these parameters will have no influence on the location of 
the iso-A:^o lines (since neither a nor Yq vary). For given physics, the knowledge of the 
effective temperature, the luminosity and the average large separation of the star helps 
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Table 3. Difference of effective temperature and luminosity between sample models, AC 
or OM, and new models with Azzq=90 /tHz but with a variation of either the initial 
helium fraction or the mixing-length parameter. The effective temperatures are given in 


Model 


Aq 

AYo 

-0.05 

-h0.05 

-0.01 

+ 0.01 

AC 

ATefJ 

-5.5 

+26.2 

-53.1 

+69.8 


A(L/L0) 

-0.027 

+0.029 

-0.079 

+0.142 

OM 

ATeff 

-7.3 

+9.7 

-57.4 

+66.3 


A(L/Lg) 

-0.009 

+0.016 

-0.073 

+0.084 


Table 4. Variation of the effective temperature, luminosity and average large separation 
for models AC and OM when subjects to a small variation of mass or age. The mass 
are given in M©, the age in Myr the effective temperatures in K and the average large 
separation in /rHz. 


Model 


A M 

At 

-0.01 

+ 0.01 

-100 

+ 100 


ATeff 

+3.3 

-1.5 

+1.2 

+1.8 

AC 

A(L/L0) 

+0.003 

-0.003 

+0.007 

-0.006 


A(A!/o) 

-3.0 

+3.0 

+1.0 

-1.0 


ATeff 

+8.0 

-3.1 

+11.1 

-11.1 

OM 

A(L/L0) 

-0.017 

+0.017 

-0.002 

+0.002 


A(Ai/o) 

+3.0 

-2.0 

+1.0 

-1.0 


to constrain values for its mass and its age. To estimate the precision obtained on both 
of these parameters, we inversely computed models including the same physics, but in 
which we changed either the mass or the age. Table01shows the effect of a mass variation 
of ± 0.01 Mq and an age variation of ± 100 Myr (only one parameter varying at once) 
on Teff, L/Lq and Ar'o- 

It appears that in all cases it does not lead to major variations, neither in effective 
temperature nor in luminosity. In contrast, variations in Avq are much more significant 
and do not agree with the seismical constraints in case of a mass change of the order of 1%. 
It thus gives an upper limit of 0.01 Mq for the precision obtained with asteroseismology 
on the mass of the star. 

In case of an age shift of ± 100 Myr, the resulting models have average large sep¬ 
arations corresponding to the limit values we used. Consequently, we expect, for given 
physics, an accuracy of the order of 100 Myr on the age of the star. 
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5. Discussion and conclusions. 

For the first time stellar oscillations have been detected and an asteroseismic analysis has 
been performed for a planet-hosting star. With a magnitude V=5.1, ^ Ara also represents 
the faintest solar-like star with detected p-modes. The observing run was a real success, 
as 43 modes could be identified, leading to a clear echelle diagram with large separations 
of 90 /rHz. Note that we also discovered during this run the smallest exoplanet ever 
observed, with a mass of 14 earth masses. 

The seismic analysis of this star leads to precise constraints on its position in the 
HR diagram, inside the observational box obtained from spectroscopy and the Hipparcos 
parallax. 

The ultimate aim of this run was to try to determine whether the star is completely 
overmetallic or only in its outer layers. In the first case it would suggest that the star 
and its planets formed out of an interstellar nebula which was already metal rich, while 
in the second case it would mean that accretion at the beginning of planetary formation 
was much more active than generally thought. 

We found an interesting test to decide between these two fundamental scenarios : the 
overmetallic models which can account for the observable parameters of the star have a 
convective core while the accretion models do not. We chose two of these models to go 
further in the analysis, by comparing the large and the small separations. We found that 
the presence of the convective core in the overmetallic model presented a signature in 
its deviations to the asymptotic theory : the lines I = 0 and I = 2 cross in the echelle 
diagram, which does not happen in the accretion model. 

Unfortunately this signature clearly appears around frequencies of 2.5 mHz which 
correspond to the highest frequency we were able to detect. At first sight, the graphs 
which display the observed and computed small separations (Figures 0 and 0) seem in 
favor of the accretion scenario, but we definitely need more observations to go further. 
It would be necessary to increase the signal to noise ratio for frequencies around 2.5 to 
2.7 mHz, which lie in the crucial region where the small separations for the overmetallic 
and the accretion models clearly diverge. 

In spite of a large amount of observational constraints on the bright star fi Ara, we 
were not able to find a unique solution for the mass of the star and propose two possible 
masses which differ by more than 10%. Further additional Doppler observations, in mono 
or multi-sites, will not easily resolve the question. In the frequency range we are trying to 
reach, p-mode amplitudes are lower than 10 cm.s“^ and we can suspect that, due to their 
short life time, it will be very difficult to extract them from granulation noise. Future 
interferometric determinations of the radius of /r Ara could help choosing one of the two 
scenarios, in case it can reach an accuracy better than 1%. Our analysis clearly illustrates 
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the difficulties of characterizing the internal structure of stars even with a large number 
of identified p-modes. 
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